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Abstract 

The stress-energy tensor of the quantum vacuum is studied for the 
particular case of quantum electrodynamics (QED), that is a fictitu- 
ous universe where only the electromagnetic and the electron-positron 
fields exist. The integrals involved are ultraviolet divergent but it is 
suggested that a natural cut-off may exist. It is shown that, in spite 
of the fact that the stress-energy tensor of the electromagnetic field 
alone is traceless (i.e the pressure P equals 1/3 the energy density 
u), the total QED tensor is proportional to the metric tensor to a 
good approximation (i. e. P = -u). It is proposed that there is a 
cosmological constant in Einstein equation that exactly balances the 
stress-energy of the vacuum. It is shown that vacuum fluctuations 
give rise to a modified spacetime metric able to explain dark energy. 
Particular excitations of the vacuum are studied that might explain 
dark matter. 


1 Introduction 

1.1 Quantum vacuum, cosmological constant and dark 
energy 

Quantum field theory predicts that the vacuum is not empty, but filled with 
interacting quantum fields. A calculation of the energy of the fields gives 
divergent quantities, positive for free Bose fields and negative for free Fermi 
fields, as will be illustrated in section 2 below for the particular case of 
quantum electrodynamics (QED). Nevertheless the energy of the vacuum 
is irrelevant in most quantum calculations not including gravity. In fact it 
may be ignored defining the zero of energy at the level of the vacuum. In 
practice this is achieved using the “normal ordering rule” for the creation 
and annihilation field operators. 

However a big difficulty appears if we want to take the vacuum gravity 
into account, because in this case the zero of energy cannot be fixed at 
will. Two alternative solutions have been proposed, none of them completely 
satisfactory: 1) Assuming that the vacuum fields represent a mathematical 
artefact of the quantum formalism and only the excitations of the fields above 
the vacuum contribute actual energy, 2) The vacuum fields have energy, 
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but positive and negative contributions may cancel to each other, at least 
approximately. In this paper the first possibility is rejected and the second 
one is studied searching for clues for the solution of a number of well known 
difficulties. 

If there is an actual energy density of the vacuum there should be also 
pressure terms, that is a full stress-energy tensor. Indeed if the spacetime 
background is Minkowski, at least approximately, then it seems that the 
vacuum properties should be Lorentz invariant whence the vacuum stress- 
energy tensor would be proportional to the metric tensor. As a consequence 
the effect of that tensor would be equivalent to a cosmological term in the 
Einstein equation of general relativity. The existence of such a term might 
explain the observations in cosmology, which have shown that the universe 
is in accelerated expansion[1] [2], Actually the cause of that expansion is 
unknown and it is usually named “dark energy”, but it is common wisdom to 
interpret dark energy as an effect of the quantum vacuum. Other alternatives 
have been proposed that will not be commented here. 

The hypothesis that dark energy is an effect of the quantum vacuum 
poses a well known problem[4], In fact it is plausible that the vacuum energy 
density, p DE , should correspond to a combination of the universal constants 
c, It, G, that is 


Pde ~ £2^ ~ lo97k 9/m 3 . 


( 1 ) 


This is Planck’s density whose value is about 123 orders greater than the 
required dark energy density, with current empirical value [3] 


Pde = ~Pde - (6.0 ± 0.2) x 10 ~ 27 kg/m 3 , 


( 2 ) 


Pde < 0 being the pressure. The first eq.(j2J) is consistent with the vacuum 
stress-energy tensor being Lorentz invariant. 

1.2 The relevance of the vacuum fluctuations 

In this paper a solution is proposed to the problem of the big disagreement 
between the “theoretical” value eq.dH), and the empirical result eq. ([2D - I 11 
order to understand the logic of the proposal I start pointing out the relevance 
of the fluctuations of the vacuum fields, that I believe has not been fully 
appreciated in the analysis of the gravitational effects of the vacuum. For 
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the sake of clarity I will begin with arguments involving classical (rather than 
quantum) fields and Newtonian gravity (rather than general relativity). 

Let us assume that the vacuum consists of a stationary, homogeneous 
and isotropic set of interacting fields, some of the fields contributing positive 
energy and other fields negative energy and similarly for the interactions. If 
we want to escape from the absurd assumption that the field energy density 
is at the Planck scale (i. e. eq.(pQ)) a plausible hypothesis is that positive and 
negative contributions cancel to each other. A partial cancelation giving the 
result eq. ()2|) , that is more than hundred orders smaller than Planck’s density, 
looks conspiratory[4]. Therefore it is plausible that the cancelation is exact, 
which I will assume at this moment (this assumption will be later modified, 
taking quantum theory and general relativity into account, see section 4 
below). If this is the case it seems that dark energy remains unexplained. 
However I will argue in the following that, even if the average energy density 
of the vacuum is zero, the fluctuations might explain dark energy. 

In fact let us suppose that the vacuum fields fluctuate. That is, at a given 
time there are regions with energy density above the mean (that is positive 
fluctuations) and other regions with density below it (negative fluctuations). 
Thus, assuming that the mean vacuum energy is zero, the fluctuation of the 
density would be a function, p (r ,t ), Milling the condition that its space 
average is zero, that is 


lim / p (r,t) d 3 r = 0. (3) 

Jv 

We must assume, as in special relativity, that energy gravitates. Then we 
should accept that Newtonian gravity associates negative (positive) gravita¬ 
tional potential to positive (negative) mass-energy. (Of course Newtonian 
gravity is not consistent with special relativity, but for a simplified argu¬ 
ment I may combine both). Thus, using Newtonian theory, the gravitational 
energy, Eq, of the vacuum fluctuations in the volume V would be 

Eg = 7; J p(ri,i) $(ri,f)d 3 ri, (4) 

*(ri,t) = —G f P(l ~ 2,t) d 3 r 2 , 

JV l r 2 - r l| 

where G is Newton’s constant and V a volume large in comparison with 
the range of the fluctuations. As the fluctuating density p has zero mean, 
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the gravitational potential of the fluctuations <f>, which is linearly related to 
p, is also nil on the average. Thus the vacuum fluctuations would produce 
no average gravitational held at any spacetime point. This might suggest 
that vacuum fluctuations would not lead to spacetime curvature when we 
pass from Newton's to Einstein's gravity. If this is the case the density 
fluctuations could not explain dark energy. However the suggestion is hawed 
due to the fact that the equations of Einstein’s theory are nonlinear, in 
contrast with Newton's, a point that will be explained in more detail in 
section 4 below. 

There is another argument supporting that vacuum huctuations, with 
zero mean, may produce a long range effect. In fact, the gravitational energy 
E g due to the huctuations, eq. (141) . is not linear but quadratic in the fluctuat¬ 
ing density. Therefore the mean gravitational energy density, p G , needs not 
be zero. 

The quantity p G is defined by 


r e g 

Pg = J im -rr 
V —>oo V 

= — lim 

V—too 


wi pMd3r \ 


l r 2 


1-1 


p (r 2 ,t) d 3 r 2 


'v 


Changing from the variables {iq, r 2 } to the new ones {r 0 , r} via 


ro = 


r i + r 2 


", r = r 2 — ri, 


we get 


Pg = 


G 

2 


-C ( r ) d 3 r = 
r 


2v tG 


C (r) rdr, 


C (r) = 


lim 

V—>-oo 


— / p (r 0 +r/2, t) p (r 0 — r/2, t) d 3 r 0 


>v 


(5) 


We see that the mean graviational energy of the huctuations depends on the 
two-point energy density correlation, C (r), and it may be different from zero. 
The two-point correlation at equal times should depend only on the distance 
r = |r 2 — iql , as is consistent with the assumed homogeneity and isotropy of 
the vacuum. Also it fulfils the condition 


C (r) d 3 r = lim 

V—voo 


= lim 

V—voo 


1 

V 
1 

V 


/ d 3 r p (r 0 +r/2, t) p (r 0 — r/2, t) d 3 r 0 
'v Jv 


p(ri,f) d 3 ri x / p (r 2 ,f) d 3 r 2 


'v 


'V 


= 0 , 
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taking eq.flS]) into account. 

The conclusion from our classical analysis is that, in addition to the mean 
energy density of the vacuum, if any, there should be a gravitational energy 
density caused by the density fluctuations. In section 3 I will calculate the 
two-point correlation in quantum electrodynamics. In section 4 the relation 
of this correlation with the dark energy density p DE , eq.(J2]) , will be studied. 
Indeed we shall see that p G may be taken as an estimate for p DE . 

The suggestion that quantum vacuum fluctuations may give rise to an 
effective cosmological constant was made by Zcldovich[5] in 1976. The pos¬ 
sibility that these fluctuations are at the origin of the dark energy has been 
explored recently[6] using a simplified model involving general relativity. 

2 Stress-energy tensor of the vacuum fields 

2.1 Energy density and pressure. Mean and two-point 
correlations 

Our aim is to study the gravitational effects of the vacuum. According to 
general relativity the quantity relevant for that purpose is the stress-energy 
tensor. If we assume that space-time is Minkowski (which is a good approxi¬ 
mation for our present purposes, but see section 3) then stationarity, isotropy 
and homogeneity of the vacuum lead to a metric tensor defined by just two 
quantities, the energy density p and the pressure P. That is the metric tensor 
should be 

r , 1 = r| = rf = p, r 0 " = - P , r; = o if ^ v . 

If in addition we assume Lorentz invariance, then we should have P = — p. 
In this case the stress-energy tensor is proportional to the metric tensor. 
Actually I shall not assume Lorentz invariance from the start, but should 
derive it as a result of the calculations. In order to obtain the gravitational 
contribution of the vacuum it is enough to get the mean energy density, 
the mean pressure and the relevant two-point correlation functions. The 
average energy and pressure will be calculated in section 2 and some of the 
correlations in section 3. 

In the quantum context vacuum expectations of the appropriate operators 
should be substituted for averages.The expectation of the energy density is 

p — (vac |p (r, t)\ vac). (6) 
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where the operator p (r, t) is the quantized Hamiltonian density, which should 
include all quantum helds and their interactions, excluding the gravitational 
one. In this paper the gravitational held (that is the spacetime curvature) 
is considered as radically different from the remaining fields of nature, al¬ 
though it should be quantized in some form (see below, section 3.2). In the 
Heisenberg picture the Hamiltonian (or energy) density operator may depend 
on position and time, but the state should not. The expectation p, eq. (EJ) , 
cannot depend on (r, t) due to the traslational and rotational invariance of 
the vacuum. 

The two-point density correlation is the most relevant quantity for the 
study of the quantum fluctuations. It is defined in quantum field theory by 


C (ri, H; r 2 , t 2 ) = (vac\p p (r 2 ,t 2 )\vac), (7) 

provided that p (rq, t\) commutes with p (r 2 , t 2 ). If they do not commute the 
correlation so defined might not be real, having an imaginary part. If this is 
the case we will define the correlation in the form 

c = - (vac |p (ri, H) p (r 2 , t 2 ) + p (r 2 , t 2 ) p (iq, H)| vac ). (8) 

If we integrate the Hamiltonian density operator, p (r ,t) , with respect to 
r we get the Hamiltonian operator 


H — f p (r ,t) d 3 r. (9) 

Jv 

I will consider the volume V as finite (e. g. a cube of side H 1 / 3 ), but take 
the limit V becoming the whole space at the end of the calculations. I 
shall assume that the vacuum state is the eigenstate of H with the minimal 
eigenvalue, which I will label E vac . That is 

H | vac) = E vac | vac), E vac = Vp vac , 

with no other eigenvalue of H smaller than E vac . Hence the r 2 integral of 
the two-point correlation gives 


C (ri,H; r 2 ,t 2 ) d 3 r 2 = ( vac 


>v 


P ( r i, ti) / p (r 2 , t 2 ) d 3 r 2 


'v 


vac 


= {vac 


p(v 1 ,t 1 )H 


vac 


= {vac |p (ri,H) E vac \ vac) = Vp 2 vac , (10) 
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where eqs.([7D and ([HD have been taken into account. This equation contra¬ 
dicts eq.(??) , the reason being that here we have used a different definition, 
not making the assumption p vac = 0. Actually the relevant quantity is the 
two-point density correlation relative to the mean, therefore I will redefine 
the correlation after subtracting the square of the mean density, that is 

C (ri, ti; r 2 , t 2 ) = C M (ri, H; r 2 , t 2 ) ~ p 2 vac ■ (H) 

where C oid is given by eq.flsp . This correlation does fulhl eq.(??). 

A calculation of the energy density and the two-point correlation of all 
quantum fields existing in nature (possibly those of the standard model of 
high energy physics) is a formidable task, out of the scope ot this paper. Here 
the study will be restricted to quantum electrodynamics (QED) in order to 
get clues for the solution of the problems associated with the gravity of the 
vacuum fields. That is I will consider a fictitious world where the existing 
quantum fields are the electromagnetic one and the Dirac held of electrons 
and positrons but no other fields. 

2.2 The free electromagnetic field. A paradox 

We will work in Minkowski space, a good enough approximation for the 
calculation of the mean vacuum energy density. (But the approximation is 
not valid when we take the two-point correlation into account, as will be 
discussed in section 3.2). In the present subsection I will calculate the stress- 
energy tensor of the free electromagnetic held of the vacuum, ignoring all 
other vacuum helds. I shall assume homogeneity and isotropy, so that the 
stress-energy tensor is defined by just two parameters, density p and pressure 
P. However, as we shall see, it is not possible to assume Lorentz invariance, 
that is P = —p. 

The quantization of a classical held involves promoting the amplitudes of 
the plane waves expansion to become (creation or annihilation) operators. 
For instance the classical free electromagnetic held, in the Coulomb gauge, 
may be represented by an expansion in plane waves of the vector potential, 
A(r,i), that is 



x [ctk,e£exp (ik.r—ikt) + a£ £ e:exp (— ik.r+ikt)] , (12) 


A(r,t) 


where k — |k| . From now on I will use the notation of the book of Sakurai|7] 
(except otherwise stated) and natural units h — c — 1. However sometimes I 
will write explicitly Planck's constant h, for clarity. From eg. (TT2|) it is easy 
to get the electric held, E = —dA/dt, and the magnetic held, B = V x A. 
The polarization vector e depends on k (in fact k ■ e =0) and it may have 
two possible values so that we should write Sj (k) , j = 1, 2, but I will use a 
simplified notation except when some confussion might arise. 

In the quantized held an annihilation operator dk. e is substituted for the 
amplitude a^ £ and a creation operator a^ e , for the amplitude , whence 
the electric and magnetic helds, E and B, become vector operators. They 
may be written as an expansion in plane waves taking a quantized counterpart 
of eq. (fT2]) into account. From these expansions, that I do not write explicitly, 
it is trivial to obtain the energy density operator. We get 


Pem = - ( E z + B^ ) _ — ^ Ab(dfc/ T dk^dj^/ £ /) 


k.e k',e' 


4V 


EE -d^dk.edk'je' T A 2 dk £ dk/ £ i J , 


k.e k',e' 


where the functions K\ and K- 2 are c-numbers (not operators) given by 


A'i = Vkk'e * e exp [i (k — k') .r —i (k — k ') t] (13) 

K 2 = Vkk'e * s' exp [i (k + k') .r —i (k + k') t] , (14) 

For notational simplicity I have introduced the following “star product” 


£ * s' = e ■ e! + — [(k x e) • (k'xe 7 )], 
kk' 

with the properties 

* £ ") = £ i* e "i £l * £ J = 4 




^(£ t * e ') 2 = 2 




1 + 


k - k' 

kk' 


(15) 


It is convenient to write all terms of Pem ' n normal order, that is the 
annihilation (creation) operators to the right (left). Taking the commutation 
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rules of the field operators into account we obtain 


Pem 

Pemi 

PEM‘2 


Pemo + Pem i + Pemi-> 

2^7 

k,e k',e' 

^ee( /^2^k,£^k',£ / “1“ -^2^k^^k'jS 7 ^ ’ 

k.e k',£' 


(16) 


where p EM0 is a numerical constant (times de unit operator) given by 


Pemo 


1 

2V 


5 > 


The Hamiltonian is obtained by performing a space integral of the energy 
density eq. (fT6|h that is 


Hem = lim / p EM (r)d 3 r 
V^oo J v 

= ^fc(dJ ce Qi kie + ^). (17) 

ks 

We see that the integral cancels de V denominator and removes the spacetime 
dependence. 

For the free electromagnetic field the vacuum state, | 0), may be defined 
as the state with the minimal energy, that is the smallest eigenvector of the 
operator eq. flTTD . It is a state with zero photons and it has the properties 

«k, £ | 0) = 0, (0 | al e = 0. 

The state | 0) is different from the actual QED vacuum, | vac), that takes the 
interaction with the electron-positron field into account. It will be studied 
below. For the purely electromagnetic vacuum state | 0) the expectation of 
the energy density is 

(0 \Pem\ 0) = Pemo = y hk , (18) 

k,£ k 

where the latter equality derives from the two possible polarizations. We 
have taken into account that the two latter terms of eq. (HBD do not contribute 
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to the expectation value because either the annihilation operator placed on 
the right or the creation operator on the left gives zero when acting on the 
vacuum state. In the limit V —> oo we obtain 


Pem o 


1 

V 


l, J 


2tt 2 


k 6 dk = 


k (27r) 3 d 3 k, 

k 4 


8tt 2 ’ 


(19) 


where we have introduced a cut-off, k A iax , in the wavevectors (or the momenta, 
what is equivalent in our units). We see that the integral is strongly divergent 
in the limit fc max —* oo, but it is plausible that a cut-off could be originated 
by fluctuations of the spacetime metric, as will be discussed in section 3.2. 

In order to define the stress-energy tensor we need the pressure. As is well 
known the pressure is 1/3 times the energy density for any isotropic electro¬ 
magnetic radiation. Therefore the pressure of the vacuum electromagnetic 
field would be one third the quantity eq. (H9j) . Hence the stress-energy tensor 
of the free electromagnetic field of the vacuum is diagonal, and its nonzero 
components are 

k A k 4 

„ _ max d _ max 7on\ 

PEMO ~ -g^2-> ^EM 0 - 24^2 > 

That is the tensor is traceless, fulfilling 


— 3-Pemo Pemo ~ 0) 

for any choice of the cut-off wavevector k max . 

This gives rise to a paradox because it is plausible that the stress-energy 
tensor of the vacuum is proportional to the metric tensor, as a consequence of 
the Lorentz invariance that we should assume for the vacuum (in Minkowski 
space). Equivalently the tensor in mixed coordinates should be diagonal 
with a pressure equal to minus the energy density. Therefore if the vacuum 
fields are real, which is the fundamental hypothesis of this paper, there is 
a contradiction between the stress-energy tensor predicted by the Maxwell 
equations (for any free radiation field) and the Lorentz invariance of the 
vacuum fields. Actually it may be argued that neither the density nor the 
pressure are well defined, eqs. (l20j) being divergent (or cut-off dependent), but 
there is a problem in any case. I think that the plausible solution is that, as 
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the vacuum contains many fields, the stress-energy tensor of all interacting 
fields toghether is well defined and Lorentz invariant , but the tensor of every 
free field alone is not. (The divergences might be removed by an effective 
cut-off, see section 3.2 ). Studying the stress-energy tensor of all fields is 
beyond the scope of this paper, but a support to our hypothesis results from 
the study of the electron-positron field to be made in the next subsection. 


2.3 The electron-positron vacuum field 

We need the stress-energy tensor of the Dirac field, that is the energy density, 
p, and the pressure, P. In order to avoid any confussion caused by the 
different notations used in the literature, I will write p and P in terms of the 
original Dirac matrices ck*, and ft (rather than the gamma matrices that have 
been defined in several different forms, compare e. g. [7] and 0 )- Thus the 
corresponding operators may be written 


Pd ~ 


p(k_ 

r D ~ 


1 

2 


- V — 

° 1 v dt 


if a , 


dxi 


dt * J ’ 

(21) 

dif* 

dx> ’ 

(22) 


where if and if are quantized fields and j may be either 1, 2 or 3, the resulting 
pressure being the same in the three cases due to the assumed isotropy of 
the vacuum. Expanding if and if* in plane waves we get 


^(r,t) 



exp (ip ■ 


r — iEt ) 



exp (—ip ■ r + iEt ), 


^ "j" /I ^ s\ 

if (r, t) = y y z2J-fij b U u U e xp(-ip-r + iEt) 

p,s V 

T H exp (*p ' r ~ iEt ) > 

p,s V 
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where b p>s yd P ,sj is the annihilation operator of an electron (positron) with 

momentum p and spin s(= 1,2), and b^ p s the corresponding creation 

operator and u,vf,v, id are appropriate spinors. (Our notation follows the 
book of Sakurai[7].) Inserting these expressions in eq. (12T)l we get p D as a 
sum of 4 terms with a product of 2 operators each. We want all the terms 
with the creation and annihilation operators in normal order so that, taking 
the anticommutation rules into account, we should make a replacement as 
follows 


d ps exp (ip • r — iEt ) d p , s , exp (—ip'r + iE't ) 

= hpp'<5 ss ' — d]y s , exp (— ip'-r + iE't ) d ps exp (ip ■ r — iEt ), 

where h pp / and S ss ' are Kronecker’s deltas. Thus we obtain 5 terms, all in 
normal order, that is 

Pd Oh t) = p m + p b (r, t) + p d (r, t) + p bd (r, t) + p\ d (r, t) . (23) 

The latter 4 terms of eq. (123 p are as follows 


Pb (r, t) 
Pd (r, t) 
Pbd (r, t) 
Pbd ( r > 


y B ( p ’ s ’ p ’ s "> K's'Ks ex P 0 (P - p') ' r-i {E - E') t }, 

pp'ss' 

y D ( p > s > Pj d P's' d ps ex P [* (P - P') ' r-i (E - E') t ] 

pp'ss' 

^2 F (p,s, p', s') d ps b p i s t exp [i (p + p') • r-i (E + E') t ], 

pp'ss' 

y F * ( p ’ s ’ p ’ s ') b l>'s' d U exp H (P + P') ‘ r+i ( E + E ') t\ ■ 

pp'ss' 


Getting the functions B , D, F and F* is straightforward but we report only 

i 12 

the expression of | F\ after the spin sum, which will be used in section 3. We 
get 


Ei f i 


m 


4 EE' 


-(E - E') 2 ^2\vl(p)u s/ (p')| 


^ (E ~ E’f 


1 + 


p ■ p' — m 2 
~EE' 


(24) 
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where we have taken into account the expressions of the spinors u and r; [T1J. 
This leads to 


K (p) u °' (n)f 


(E + m) ( E' + m) p q 


2 m 2 E + m E' + m 


ss' 

whence putting 


p 2 = E 2 — m 2 = (E + m ) (E — m ), 
q 2 = E' 2 — m 2 = (E' + m ) (E* — m) , 


it is easy to get eq.(12D . 

The first term in eq. (!23|) . p D0 , is a c-number, not an operator (more 
properly it is proportional to the unit operator), that is 



(26) 


Integration of eq. fl23|) with respect to r gives the Hamiltonian of the free 
electron-positron field, that is 



(27) 


The vacuum state, | 0), of the free held corresponds to the eigenvector of the 
Hamiltonian eq. ([2T]) with the smallest eigenvalue. It consists of zero electrons 
and zero positrons. Therefore from now on we will define | 0) to be the 
QED unperturbed vacuum state, which is a simultaneous eigenvector of both 
Hamiltonians eqs.( lT7|) and (127)) . Thus the state | 0) is defined, from now on, 
as having zero photons, electrons and positrons. It should be distinguished 
from the physical vacuum state, | vac), which is an eigenvalue of the total 
Hamiltonian, including the interactions. 

Our aim is to obtain the energy density, that is the expectation value in 
the vacuum of the density operator eq. (1271) . Neglecting the field interactions 
(to be studied in the next subsection), we get the energy density of the un¬ 
perturbed electron-positron held as the vacuum expectation of the operator, 
that is 


(0 | p D (r, t) | 0) = (0 | p m (r, t) | 0) = p m 
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The former equality follows from the fact that the four latter terms of eq. (fTTl) 
do not contribute to the expectation because either the annihilation operator 
placed on the right or the creation operator on the left gives zero when acting 
on the vacuum state. 

The quantity p D0 , eq. ()26]) . leads to the following negative divergent energy 
density in the limit V —> oo 


Pdo ~ 


(2vr 

1 


'- 3 £ 


PP ma 


Ed 3 p = — 7 r 


-2 


fPma 


\J m? + p 2 p 2 dp 


87r 2 
1 


Pm^V™ 2 + PLx + 2 Pmax) “ ^ aT g sillll 


Pn 


47r 2 


Pmax + Pmax ? ™ 2 + o™ 


1 4 1 4 f 2 pn 

— 77 i in 


m 


m 

+ O (p m 2 „J? 8 ) 


where we have introduced an ultraviolet cut-off of the momenta, p max - I point 
out that the negative value might be anticipated by inspection of the Hamil¬ 
tonian eq. (12711 . I 11 fact the density p D0 is just the (unperturbed) vacuum 
expectation of the Hamiltonian divided by the volume V. 

The stresses may be calculated as vacuum expectations of the operators 
Pp, eq. f| 22 jl . As in the case of the energy density, eq. f)23jl , we may write the 
operator P D as a sum of terms all of them with the creation and annihilation 
operators in normal order. We get 

PH (r ,t) = P£ + PH (r, t) + ff (r, t) + (r, i) + P« (r, t ). (29) 


The latter 4 terms are similar to those of eq. (123 p but I will not give their 
explicit expressions. These 4 terms do not contribute to the mean pres¬ 
sure. The former term of eq. (l29jl does contribute and, after some algebra, it 
becomes 



where the isotropy has been taken into account in the latter equality. That is 
the stresses along three orthogonal axes are equal and every one corresponds 
to the pressure. We get 


Pdo — 


12tt 2 


Pm ax - g m2 Pmax 


b 


+ Pmax + 


3m 4 


arg sinh 


Pv 


m 
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The sum, p ZPF (for ‘zeropoint field’), of the quantities eqs. ffTUj) and ([25]) is 
the energy density of the QED unperturbed vacuum state. It is not obvious 
whether we should identify the cut-off photon momentum (the same as the 
energy in natural units), k max , with either the maximum electron momen¬ 
tum p max or the maximum electron energy, E max = ^Jm 2 + p 2 ^^. With both 
choices the leading term of p ZPF becomes 

k 4 

PZPF = PeMO + PDO = ~ ^ (^max) > (30) 

That energy density p ZPF is the quantity removed in standard quantum 
calculations by the normal ordering rule. The sum of the pressures is positive, 
namely 

PzPF = PeMO + PdO = g™ a 2 X + o ax) • (31) 

More precisely, neglecting only terms that go to zero when the cut-offs go to 
infinity, we get 


Pzpf + P zpf — 


that suggests the identification 


67T 2 


[ k 


4 

max 


Pmax^maxl 


rtT . . / Pm a x 

-- arg smh - 

47 t 2 V m 


lA — r, 3 F 

"'max rmax^max- 


(32) 


It is remarkable that the leading terms give rise to a stress-energy ten¬ 
sor fulfilling the Lorentz invariant relation P ZPF = — P ZPF , although every 
vacuum field alone does not fulfil that relation. On the other hand the sum 
eq. (1521) has positive pressure and negative energy density contrary to what 
we might expect. In section 4 I will propose a solution for this anomaly. The 
relation eq. (|32|) corresponds to the neglect of the interaction, to be studied 
in the next subsection, and it has been calculated only for the particular case 
of QED. Therefore no strong conclusion may be obtained, but our result sug¬ 
gests that the stress-energy tensor of all interacting fields, might be Lorentz 
invariant in the vacuum (in Minkowski space). 


2.4 Corrections for the interaction 

The physical vacuum of QED, | vac), is different from the zeroth order vac¬ 
uum, | 0), studied above. The latter is the eigenvector, with the smallest 
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eigenvalue, of the unperturbed Hamiltonian H 0 = H EM + H D , see eqs. (IT71) 
and (12711 . The former is the eigenvector of H = H EM + H D + H int with 
the smallest eigenvalue, that is taking the interactions into account. Find¬ 
ing | vac) as an exact eigenvector of H is not possible and we should use 
a perturbation method, that is to calculate it as an expansion in powers of 
the coupling constant, the positron charge, e. Actually only even powers of 
e would appear and the result becomes an expansion in powers of the fine 
structure constant a = e 2 /(Antic) ~ 1/137. 

The interaction Hamiltonian (or energy) operator may be written, in the 
Coulomb gauge, 

Pint (?,t) = -e-0 oit/j • A. (33) 


The operators Pp,Pj? and A contain two terms each when expanded in plane 
waves, each term corresponding to an infinite sum. One of these terms has 
creation operators and the other one annihilation operators (see e. g. eq. (TT2]l 
for the electromagnetic potential vector). This gives rise to 8 terms for p int , 
eq. (153jh Integrating with respect to r inside the volume V leads to the 
interaction Hamiltonian. Only two terms survive and we get 


Hint 


- E 


m 


V 3 / 2 V2kEE 


:U 


t 


p) « • £V S ' (q) 


p,q,k,s,s',e: 

xy +q ,koi,«!>p,. <i h' + . 


(34) 


where h p + q ,k is a Kronecker’s delta, h.c. means Hermitean conjugate and 


E = \/p 2 + m 2 , E' = \Jq 2 + m 2 . 


One of the terms of the Hamiltonian may create triples electron-positron- 
photon and the other term may annihilate triples. The integral in r causes 
that the Hamiltonian eq. ('13411 is invariant under traslations and rotations, 
whence it can couple only states with the same total momentum and total 
angular momentum as shown below, see eqs. (135j) and (|SSD • 

The physical vacuum, | vac), may be obtained as an expansion in powers 
of the coupling constant a, but for the calculation of the energy to first order 
we need | vac) only to order ct 1//2 , that is we should not consider states with 
more than three particles (i. e. one electron, one positron and one photon). 
We have to that order 


vac) = N 


°) + c i I •?) ’ 

o 


(35) 
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where | j) represents a state with one photon-electron-positron triple (differ¬ 
ent j correspond to the particles having either different momenta or spins or 
both). N is a normalization factor. The coefficients Cj may be got from the 
interaction Hamiltonian, that is 


c ; = 


Hj 


int 


En — Ej 


i E j = \J 


Ho 


j)=VPr 


(36) 


By the properties of H int , see above, the states | j) have the same net mo¬ 
mentum and angular momentum than the vacuum, that is zero. However 
they possess positive energy above the vacuum state | 0). 

The three-particles state may be got from | 0) by the action of creation 
operators as follows 

I i) = h iJi.A lS ' I °>,i = ( k A, p,s, q,«'} , (37) 


with the constraints 


k + p + q = 0,£-|-s-|-s' = 0. 

The unperturbed energy of one of these states is the corresponding eigenvalue 
of the unperturbed Hamiltonian, that is 

Ej = k + \Jm 2 + p 2 + \/m 2 + q 2 + E ZPF , (38) 

where E Z pf is the vacuum energy in the volume V, that is the density 
Pzpfi eq.(J20]) , times the volume V. Corrections of higher order in a would 
involve a sum similar to eq. (!5!)|) , but containing more general states, that 
is states having n triples electron-positron-photon, n = 1,2, 3, ...All those 
states, in addition to having nil momentum and angular momentum, possess 
the quantum numbers of the unperturbed vacuum state, | 0), i. e. zero 
electric charge (and zero leptonic number). Thus the vacuum state | vac) 
also possesses those properties. 

There are other states defined by a sum similar to eq. (ES5D , but orthogonal 
to the vacuum state | vac). That is states of the form 

I vac {l ) = | j), 

j 
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where now | j) represents a general state with n triples electron-positron- 
photon, j — 0 corresponding to the unperturbed vacuum state | 0). Thus for 
clarity I will rewrite the vacuum state eq. (I35[) in the form (not normalized) 


vac) = 


) = H c i \j)- 


Amongst the states | vac orthogonal to | vac), that is fulfilling 

V ¥i = o. 


the one with the minimal energy will be named “first excited vacuum state”, 
labelled with l = 1. Similarly the “second excited vacuum state”, l = 2, 
will be the state with minimal energy amongst those orthogonal to both the 
vacuum state and the first excited vacuum state, and so on. 

ft is common practice in quantum held theory to name “excitations of 
the vacuum” the states with particles. For instance a state with one electron 
or a state with one electron and one photon, etc. However the states here 
called “excited vacuum states” are a particular kind of vacuum excitations, 
namely those characterized by having the same quantum numbers of the 
vacuum (e. g. zero charge if we restrict ourselves to QED) and also zero net 
momentum and net angular momentum. An interesting question is whether 
these states appear somewhere in nature. My conjecture is that, with some 
slight modification, they may be the constituents of the “dark matter”, as 
will be commented in section 5. 


3 Two-point correlations of the vacuum fields 

Vacuum fluctuations may be characterized by the two-point correlation of 
the vacuum stress-energy tensor. Assuming homogeneity and isotropy that 
tensor contains at most two independen parameters, namely energy density 
and pressure. However the correlations of the different components of the 
stress-energy tensor may depend on the direction. E. g. the correlation 
(T xx (r a )T yy (r&)) may depend on the angles of the vector r a — with the 
directions X and Y. In this paper I will not calculate the two-point cor¬ 
relations of the vacuum fields in detail. But for illustrative purposes I will 
get the correlation of the density at equal times, which depends on a single 
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parameter, namely |r a — r&|. I shall do that in the approximation of zeroth 
order, that is involving quantum expectations in the unperturbed vacuum 
state | 0), rather than the physical vacuum state | vac). 

3.1 Two-point density correlation for free vacuum fields 

In this subsection I will work in Minkowski space, but see next subsection 
for a discussion about the validity of this assumption. Using Cartesian coor¬ 
dinates I shall calculate the correlation eq. flTTT) , that is 

C(r) = (vac |p (ri) p (r 2 )| vac) — p 2 , p = (vac |p| vac) . (39) 

The energy density operator is the sum of 3 terms, namely 

P = Pem + Pd + Pint j (40) 

that have been discussed in section 2, see eqs. (TT6l) and (1231) and (1331) . 

I start proving that a numerical constant added to the energy density 
operator does not change the value of C(r). In fact, if we substitute p + K 
for p in eq. (13U1) we get 

C'(r) = (vac | [p (ri) + K] [p (r 2 ) + K] \ vac) — (vac \ [p + K] \ vac) 2 

= C(r ) + K [(vac \p (rx) | vac) + (vac \p (r 2 )| vac)} — 2K (vac |p| vac) 
= C(r), 

where we have taken into account that (vac \ vac) = 1 and (vac\p (ri)\vac) 
cannot depend on r due to the traslational invariance of the vacuum. The 
result implies that we may ignore the term Pzpfi eq. l(30l) , which is equiv¬ 
alent to putting Pzpf = 0, for the calculation of the two-point correlation 
C(r). That is from now on I should use eq. fl40p ignoring the energy density 
operators p EM0 and p D0 

The calculation for free fields, that is to zeroth order in the interaction, 
consists of substituting | 0) for | vac) in eq.([28j) and neglecting the term p int 
in eq. f!40p . Thus the two-point correlation for free fields (i.e. the unperturbed 
correlation) is the following 

Co (r) = (0 \[Pem ( r i) + pp> (id)] [p EM (r 2 ) + p D (r 2 )]| 0) (41) 

= (0 | Pem ( r i) Pem ( r 2 ) I 0) + (0 | p D (q) p D (r 2 ) | 0), 
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where the equality is a consequence of the fact that the cross terms (0 | p EM (iq) p D (r 2 ) | 0) 
and (0 \p D (iq) p EM (r 2 )| 0) give no contribution. The electromagnetic free 
held contribution becomes 


C E Mo( r ) — (0 I PeM ( r l) PEM ( r 2 ) I 0) — (0 I PeM2 ( r l) PEM2 ( r 2 ) I 0) 
1 
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where we have taken into account that the term involving p EM1 ( r ) does not 
contribute (it annihilates the unperturbed vacuum state | 0), see eq. (??)). 
Hence putting the operators in normal order, using the commutation rules, 
and introducing an exponential convergence factor we get after some algebra 


CWo(r) = kk> 


k,e k'.e' 


1 + 


k - k' 

kk! 


-i 2 


exp [i (k + k') .y—e (k + k ')], 


where eqs. (fl4|) and f|T5]) have been taken into account. 
Substituting integrals for the sums this leads to 

C EM0 (r) = ————g [ kd 3 k [ k'd 3 k' 

4(27 t) J J 

x exp [i (k + k') .r—e ( k + k')]. 


1 + 


k - k' 

'W 


(42) 


The exponential convergence factor exp [—£ ( k + k')\ is so chosen in order 
to get easy integrals, analytical in some cases. It is plausible that the re¬ 
sult would not depend dramatically on the type of cut-off used, taking into 
acount that the approximations involved in our calculation will allow only 
obtaining a rough approximation. The details of the integals may be seen in 
the Appendix and the result is 


Cemo { r ) 


3r 4 - lOrV + 3e 4 

7r 4 (e 2 + r 2 ) 6 


(43) 


The correlation Cemo ( r ) is positive for small r as it should, because for r = 0 
the correlation becomes the variance. That is a fluctuation above (below) 
the average energy density is most probably close to another one also above 
(below) it. For r £ (^/3, 3e) the correlation is negative and it is again positive 
for r > 3 £. 
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A calculation of the two-point correlation of the free electron-positron 
field starts from eq. (1231) , whence we get for t = 0 

C D o{r) = (0 | p D (rj) p D (r 2 ) | 0) = ^0 | p bd (iq) pl d (r 2 ) | 0^ 

= ^2 1-^ (P> s 5 P’ s 0| 2 exp [i (p + pO • r], 

pp'ss' 


where we have taken into account that the terms p b (r) and p d (r) do not con¬ 
tribute because they annihilate the state | 0). Hence inserting the expression 
of |F| 2 eq. (l24|) . we obtain 


C D o(r) = 


4V 2 


Z IE ~ E 


l \2 


PP 
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p • q — m 




d A p / d 6 q (E - E 


EE’ 
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4(2vr) 6 ^ 
x exp [i (p + q) • r—e (p + q)} , 


1 + 


exp [i (p + p ) ■ rj 


p • q — m 


EE' 


(44) 


after introducing a cut-off in the particles momenta. The integrals are in¬ 
volved and they will not be calculated exactly. If m « e _1 we might 
approximate E = \/p 2 + m 2 ~ p, E' ~ q , whence eq. (l44|) leads to 


C D o (r) - 


3 r 4 — 10r 2 e 2 + 3e 4 
27r 4 (r 2 + e 2 ) 6 


+ O (m 2 ) . 


(45) 


Details of the calculation may be seen in the Appendix. 

The contributions of order a to the correlation are straightforward al¬ 
though lengthy. They will not be reported in this paper. 


3.2 The relevance of metric fluctuations: A natural 
cut-off 

If spacetime is Minkowski then the two-point energy density correlation of 
the vacuum fields may depend only on the interval between points that is, 
with an obvious metric, 


(vac | p (ri, H) p (r 2 , t 2 ) \ vac) = C (a), a 2 = (t 2 - hf - |ri - r 2 | 2 . (46) 
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This would be a consequence of the Lorentz invariance of the vacuum. Of 
course the correlation might depend on whether the interval is timelike or 
spacelike (in other words if a as defined above is real or imaginary). However 
a Minkowski spacetime is not compatible with the existence of fluctuations. 
Therefore eq. (14611 is at most an approximation. 

The proof of incompatibility is trivial and it follows. Let us consider a 
Minkowski space and two arbitrary points, A and B , with coordinates (r^, tjf) 
and (r bAb), respectively. It is always possible to find another point C with 
coordinates (re, tc ) which is lightlike separated from A and also from B. In 
fact any point in the intersection of the light cones of A and B fulfils that 
condition. Then eq. (H(l implies that the vacuum energy density in C will 
be the same as in A and the same as in B. But the points A and B being 
arbitrary we conclude that the density will be the same in all points of the 
Minkowski space. There would be no fluctuations at all! 

We may arrive at the same conclusion taking into account that if the 
stress-energy tensor of the vacuum fields fluctuates, then the metric tensor 
should also fluctuate by Einstein equation of general relativity. The argu¬ 
ment is clear in the context of classical physics and it should be taken into 
account also in the quantum realm. In the calculation of the two-point 
correlation eq. ffXTj) we have used (implicitly) the Minkowski classical metric 
diag (1, 1 , 1 , — 1) , but we should use a quantum metric tensor operator g lw (x), 
where x represents the four coordinates of a point. That operator must be 
related to the density operator, p(x), (more generally to the stress-energy 
tensor operator) via a quantized Einstein equation. However a satisfactory 
quantum gravity theory is not yet available. Therefore a correct calculation 
of the two-point correlation cannot be made. With our present knowledge 
it is only possible to perform approximate calculations that mix classical 
general relativity with quantum theory in some form. 

The standard method to deal with problems involving both quantum 
theory and general relativity is the use of the “semiclassical approximation”. 
It consists of calculating quantum averages in a given spacetime (usually 
Minkowski) and putting the averaged quantities on the right side of Einstein 
equation, this treated as classical, in order to get the spacetime metric. How¬ 
ever the semiclassical approximation is not good enough for our purposes in 
this paper. In fact its use leads to lossing two relevant properties of the 
two-point density correlations of the vacuum fields. Firstly it does not al¬ 
low deriving the existence of a natural cut-off in the energies of the particles 
associated to the vacuum fields, as is made in the following. Secondly the 
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approximation prevents calculating the long range effect on the metric of the 
vacuum fluctuations, as will be discussed in section 4. 

In order to go beyond the semiclassical approximation I will use a different 
approach that consists of treating the metric fluctuations, actually deriving 
from the quantum character of the metric, as if they were associated to a 
classical statistical ensemble of metrics. That is I will replace the correct, 
fully quantum, vacuum expectation by a classical average over a statistical 
distribution of quantum vacuum expectations, each expectation calculated 
with a given classical metric. 

More explicitly I will consider a probability distribution on an ensemble, 
J, of metrics, all of them close (in some sense to be specified later) to a ref¬ 
erence Mikowski metric diag( 1,1,1,—1). I will use for every metric of the 
ensemble J the same (Cartesian) coordinate system as in Minkowski space, 
but obviously a different metric tensor, g^ u , and a different volume element. 
Our problem is to find a more accurate substitute for the calculation made 
in the previous subsection, i. e. C 0 eq. flTTf . Specifically the calculation in 
Minkowski space required integrals involving dynamical variables, like mo¬ 
menta of the particles, multiplied by position vectors, see for instance eqs. (I42p 
or (1441) . We should find substitutes for those integrals in any metric gf u , and 
also we should find the probability distribution in the ensemble of metrics. 
The calculation would be involved, or impossible, and I shall make a simpli¬ 
fication. I will just consider a family of extremely simple metrics in order to 
get hints about the main consequences of the metric fluctutations. 

Let us consider a family of metrics with the form diag( 7 ,1,1,—1), 7 
being a (positive) real number close to unity, and let us study the change 
produced by that metric in a typical integral appearing in eq. (l42|) (where now 
we remove the cut-off, exp [— ek ],because deriving the existence of a cut-off 
is just the purpose of the present calculation ). For instance we will study 
the change produced in the following integral 



where k is the momentum (or wavevector) of a virtual photon of the vacuum. 
Actually the new metric, diag( 7 ,1,1, —1), corresponds to a Minkowski space 
too, therefore it is plausible that a calculation in the new metric would lead 
to the integral 
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where the overall factor 7 takes the change of the volume element into ac¬ 
count. Now let us assume that the distribution of metrics in the family 
corresponds to a Gaussian probability distribution, /(y), where 7 is now 
taken as a random variable. That is 

/(7)<iT= 27?; exp (- i2 4^) <i7 ' 


(the probability /(y) is only approximately normalized because the variable 7 
cannot have negative values). The average of the integrals / 7 in the ensemble 
of metrics would be 


0 

where with negligible error we have extended the 7 integration down to — 00 
and approximated 7 by unity in the overall factor. 

We could extend the family including other metrics like cliag( 1, 7 ,1, —1) 
or diag(l, 1, 7 , — 1). The net result would be that the average value of the 
integral (/ 7 ) should exhibit a cut-off in the photon momenta. This illustrative 
calculation supports the intuitive idea that randomness implies fuzzyness 
and fuzzyness causes some average over photon momenta in the exponential 
exp (ik • r) . The fuzziness increases with the momentum k, which amounts to 
an effective ultraviolet momenta cut-off. However neither the calculation nor 
the intuition allow us finding the form of the cut-off. Actually for the purpose 
of getting just a rough estimate of the effect of the vacuum fluctuations the 
form of the cut-off is irrelevant. However an estimate of the maximum allowed 
particle’s momenta (say, an estimate of the parameter cr) is required. Finding 
such an esmate is difficult, but arguments for its plausible order of mangitude 
follow. 

It is common wisdom that there are quantum metric fluctuations at the 
Planck scale. If these are the only relevant metric fluctuations then the mo¬ 
mentum cut-off should be of order the inverse of the Planck length, which 
would imply an extremely large value for the two-point correlations (e.g those 
calculated in eqs. (l42|) and (jUj) ). That is the energy of the vacuum fluctu¬ 
ations would be of order the vacuum energy itself. It seems more plausible 
that the cut-off momentum is far smaller, of order the typical particle scales. 
That is it should be of order the inverse of a Compton wavelength, h/mc , m 
being a typical (or average) particles mass. 



A/(7)^7 - 
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4 


Vacuum fluctuations and dark energy 


4.1 The need of a cosmological constant in the Einstein 
equation 

It might be that the energy density of the different vacuum fields cancel to 
each other, and similarly for the pressure, but our calculations within QED, 
sections 2.2 and 2.3, do not support that assumption. A weaker assumption, 
consistent with the desired properties of the vacuum, is that the whole set of 
vacuum fields have a stress-energy tensor proportional to the metric tensor. 
Furthermore our QED result suggests that this may be the case, but the 
proportionality constant is divergent, or huge if we assume a natural cut¬ 
off at the Planck scale. However this gives rise to the problem commented 
in the introduction section, namely a big discrepancy between theory and 
observations. A plausible solution to the problem is proposed in the following. 

We return to the discussion of section 1.2 about the relevance of the 
vacuum fluctuations, but now in the realm of general relativity rather than 
Newtonian gravity, although still within classical (not quantized) theories. 
Einstein equation may be written 

- \%.R - Ag„>. = -8 nG (r™“ + r;r), (47) 

(here I use the notation of the book by Weinberg[9j.) The latter term of the 
left side involves the cosmological constant A. On the right side the former 
term is due to matter and the latter to the vacuum fields, both stress-energy 
tensors, T™ att and T^“ c , depending on spacetime coordinates. The variation 
of the vacuum tensor takes into account that vacuum fields fluctuate, see 
section 1-2. We will assume that these fluctuations are short ranged, as 
discussed at the end of section 3.2. 

Our proposed solution is that either the mean vacuum stress-energy tensor 
is zero or it is precisely balanced by a cosmological constant, whilst the 
fluctuations give rise to a long range effect able to explain dark energy. If we 
pass from classical to quantum theory, then vacuum expectations should be 
substituted for the averages over the volume V. Thus the main assumptions 
in this paper may be stated as follows: 

Proposition 1 The quantum expectation value of the stress-energy tensor 
operator due to the vacuum fields is proportional to the expectation of the 
metric tensor (or zero). 
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That is 


(T;r) = (vac I f;“ c I vac) = T vac (vac \ g, u \ vac), (48) 

where T vac is a constant and g^ is the quantized metric tensor operator. 
This implies that (the vacuum expectation of) the stress-energy tensor of 
the combined vacuum fields has the form of a cosmological constant in the 
Einstein equation, a hypothesis that is common wisdom. That wisdom how¬ 
ever suffers from a big difficulty, namely the said cosmological constant would 
be huge, which leads us again to the problem discussed in the introduction 
section, that is the disagreement between eqs. © and d2J) • For me the most 
plausible solution to the problem is the following assumption: 

Proposition 2 There is a cosmological term in the Einstein equation that 
exactly balances the effect of the mean vacuum stress-energy tensor. That is 


8nGT vac = A. (49) 

Assuming that the cancelation is exact we avoid any claim of conspiracy 
(see section 1.2). Eq. (1491) suggests writing the quantized Einstein eq. (l47l) in 
the form 

G„„ = -SttG (f;r“ + «■;“) . 5T£° = f“° - r“g„„, (50) 

where G^ u , T™ att and 5T™ C are tensor operators. The meaning of this equa¬ 
tion is not obvious because there is not yet a satisfactory quantum gravity 
theory and therefore it is not known how to relate the metric tensor operator 
g^ v with the Einstein tensor operator G^ u . 1 shall present in section 4.3 a 
procedure to give a sense to some approximation of eq. flSO]) and to get sen¬ 
sible solutions. 1 will point out that this equation departs from standard 
semiclassical approach to general relativity, where the source term in a given 
quantum state should be the expectation value of the relevant field operators 
in that state. In fact, the semiclassical alternative to eci. dSO]) would be 

= -87rG(vac | T™ att + 5f™ c \ vac), (51) 

but the term (vac \ 5T™ C \ vac) has been assumed nil, see eq. fl48l) . Con¬ 
sequently the effect of the vacuum fluctuations would be lost in a standard 
semiclassical treatment. 
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4.2 The gravity of vacuum fluctuations. A toy model 

In order to study the spacetime curvature produced by density and pressure 
fluctuations we must start from Einstein eq. (l47l) . which leads to the quan¬ 
tized Einstein eq. flSOp with our assumptions. Ignoring the effect of matter, it 
simplifies to 


G, u = -8 ttG 5T™ c , (52) 

where the Einstein tensor operator G is a functional of the metric tensor 
operator g ltl/ . The functional G, w [^a<t] would contain g\ a and its first and 
second derivatives with respect to the coordinates, but we do not know what 
is the exact relation, something that should be derived from a quantum 
gravity theory. As a consequence eq. d^jl is not well defined. Nevertheless 
for the study of vacuum fluctuations an approximate quantum equation may 
be substituted for eq. (1521) , and sensible solutions found, as shown in the next 
section. 

Before presenting that solution, I will work an extremely simplified model 
where the two-point density correlations depend on a single parameter, the 
radial coordinate. Of course the model is far from realistic, but it provides 
hints about the actual consequences of the vacuum fluctuations. I consider a 
static problem with spherical symmetry in standard (curvature) coordinates 
{r, 9, 0, t} with some stress-energy tensor whose components fulfil 


vac {8f va % vac} = 0, 

as a consequence of eq. (1481) . Using the metric 

da 2 = A (r) dr 2 + r 2 [d9 2 + sin 2 Odcj) 2 ] — B (r) dt 2 , 
one of the components of the classical Einstein’s equation is [9] 

i(j) = 1 - SnGr2p{r) ’ 


(53) 


(54) 


(55) 


with the initial condition A (0) = 1. We shall work in quantized general 
relativity. Although no satisfactory quantum gravity theory exists yet, it is 
plausible that the relations between the quantum operators are the same as 
the classical ones whenever there is no problem with the non-commutativity 




of the operators. Therefore the quantized counterpart of eg . (1551) relevant for 
our problem would be 


(l) = 1 ~ 87rGr2 ^ ( r ) ’ ( 56 ) 


where 87T dp is the component of the stress-energy tensor, (dT vac ) v ^ derived 
from the vacuum fluctuations. I stress that the average tensor due to the 
vacuum fields does not appear in eq.(??) because 1 have assumed that it is 
exactly canceled by a cosmological constant. However for notational simplic¬ 
ity 1 will substitute p for dp from now on, but take into account that we will 
have 

(vac |p| vac ) = 0. (57) 

The exact solution of this operator equation is trivial and it may be written 
in closed form, that is 


A 



8ttG 


x 2 p ( x ) dx 


-l 


(58) 


Now we assume that Newton constant G is small (e. g. in comparison with 
ch/m 2 , m being a typical particle mass). Thus we expand eq. (l58li in powers 
of G giving 

, , 8v tG f r J 

A = 1-|-/ x p (x) dx 

r Jo 

647 t 2 g 2 r r 

+ — ~ 2 — J x 2 p (x) dx J y 2 p (y) dy + O (G 3 ) . (59) 

The quantized densities at different points may not commute, that is we 
might have 

[p 0) , P ( 2 /)] 0, 

but this is irrelevant in eq. (1591) because the latter term is itself the square of an 
operator. The interesting quantity is the vacuum expectation of the metric, 
which would provide the observable space-time curvature. The expectation 
of the above equation leads to 


vac 


A 


vac ) ~ 1 + 


64tt 2 G 2 


x 2 dx / y 2 dy (vac\p (x) p (y)\vac ), (60) 
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where eg. (15311 has been taken into account. As a conclusion we have obtained 
the expectation value of one of the metric components in terms of the two- 
point correlation of the energy density, and shown that the metric component 
may have long range (of order r) deviations from Minkowski metric, in spite 
of the vacuum expectation of the density being zero , see eq. ([57]) . I stress that 
this result could not be derived from the standard “semiclassical approxima¬ 
tion” to general relativity. In that approximation we should write, instead 
eq. (I56j) , the following 


d 

dr 


vac 


A 


vac 


= 1 — 8nGr 2 (vac \5p (r) | vac) = 1, 


whence we would get 


vac 


A 


vac) = 1 , 


that is the Minkowski value. 

We may search for an effective (classical) density, p e fj (r), providing the 
same spacetime curvature as the vacuum expectation of the metric compo¬ 
nent, eq. (160]) . It would fulfil 


vac 


A 


vac) — 1 + 


87 tG 


x 2 Peff ( x ) dx + O (G 2 ) , 


that after inserting it in eq.dfiOl) and calculating the derivative leads to 


87 tG d 

M - — Tr 


- / x 2 dx / y 2 dy (vac \p (x) p (y)\ vac) 

r Jo Jo 


( 61 ) 


This result emphasizes the main consequence to be derived from our toy 
model, namely that the vacuum fluctuations produce a metric equivalent to 
the metric derived from some effective (fictitious) density. 

The second component of the classical Einstein equation cannot be quan¬ 
tized as easily as eq. (I55|) , so that we may follow a different route in order to 
get the quantized component of the metric, B{r), see eq. (l54|) . We start from 
the well known classical solution 


B(r ) = exp 



2 M ( x ) + 87 tx 3 P ( x ) 
x - 2 GM(x) 


(62) 
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where 


M ( x ) = / 4v iy 2 p (y) dy , 

Jo 

P being the pressure (more correctly the component 5Tj of the stress-energy 
tensor eq. (|52]) divided by — 87 t). Expanding eq. (l62j) to second order in G we 
get 


B(r) = l + G 
+2 G 2 


2Mix) + 8ttx 3 P (x) , 1„ 2 

-—- —dx + -G 

x 2 


2 M[x) + 87 tx 3 P (x) 


dx 


LJo 


x 


[2 M(x) + 87 tx 3 P (x)]M(x) 


x* 


dx + O (G 3 ) . 


In order to quantize this expression we should substitute operators for the 
dynamical variables. However there is a difficulty due to the possible non¬ 
commutativity of P (x) and M(x). I propose the following rule: 

For the quantization of the product of two components of the stress-energy 
tensor we should use the symetrized product. 

This leads to 


P (x) M(x) = 




Jo 

1 r 


4t ry 2 dyP (x) p(y ) 

2t ry 2 dy P (x) p(y) + p(y)P (x) 
P (x) M(x) + p(y)M (x) 


Then the vacuum expectation of the quantized metric component B(r) be¬ 
comes 


vac 
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vac) = 1 + 4 G 2 
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vac 


M(x) 
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vac 


dx 


x* 


(63) 
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P (x) M (x) + M (x) P (x) 
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We may introduce an effective pressure, P e ff ( r ), such that 

2 M eff (x) + 8nx 3 P eff (x) 


vac 


B 


vac ) = 1 + exp 

nr 

= 87 tG 


G 


x — 2 GM e ff(x) 


dx 


x 


P e ff (x) H- In f — j p e f f (x) x 2 dx + 0(G~\64) 


It is straightforward to get P e ff(r) equating this with eq. (l63|) , but I will not 
write the result. 


4.3 Approximate solution in four dimensions 

In order to devise a method for the solution of eq. (!5Dl) , ignoring the effect 
of matter, I propose the following generalization of the method used in the 
previous toy model. 

We should write the metric tensor, gjk , as an expansion in powers of the 
Newton constant, G , that may be taken as a small parameter whence the 
deviation from the Minkowski metric is assumed small. Therefore we write 


_ (0 . n (1 . (2 
9jk - 9jk + Gdjk + G 9jk— 

Thus the zeroth order approximation of eq.(??) becomes 


G 


flU 



0, 


whence g^° k is the metric tensor of Minkowski space. The first order approx¬ 
imation will be 



where the functional derivative should be taken at g ^ = g^ k . The point 
is that this equation is linear in the unknown, therefore relatively easy to 
solve. The second order approximation would involve the second functional 
derivative, giving another linear equation, and so on. In summary we reduce 
the Einstein nonlinear equation to an infinite set of linear ones. Of course 
this is only valid because we are assuming here that the metric is close to 
Minkowski, with small fluctuations superimposed on it. As explained below 
the second order is sufficient for the study of the vacuum fluctuations. 

An advantage of the approximate method to solve Einstein eq.(??) is 
that it may be used, with plausible assumptions, for a treatment in quantized 
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general relativity. In the quantized theory the tensors T™ c and gjk become 
operators, but their relation cannot be obtained from the classical (Einstein) 
equation due to the possible non-commutativity. However, the approximate 
relations between the metric and the stress-energy tensors, involved in the 
calculation up to second order in G, are either linear or they contain products 
of no more than two components of the stress-energy operator. In the linear 
relations there is no problem of commutativity and in relations involving 
products of two operators it is plausible to use the symmetrized product. 
That is if T™ c (. x ) T^“ c (x') is a classical product, the quantum counterpart 
would be 


1/2 


rnvac 

1 (1U 


X 


rpvac 

Act 


X 


I rpvac 
1 Act 


M T. 


vac 


X) 


whee x represents the four coordinates of a spacetime point. 

The conclusions of this section are the following: 

1) It is possible to quantize the Einstein equation of general relativity in 
the post-Newtonian approximation, that is neglecting terms of order higher 
than G 2 provided we assume that the quantum counterpart of the product 
of two components of the stress-energy tensor gives a symmetrized product 
of operators. 

2) The vacuum expectation of the metric induced by the quantum vacuum 
fluctuations is the same as the one due to a classical effective stress-energy 
tensor, that is an effective mass density, p e ff, and pressure, P e ff if the space- 
time is Minkowski (at least with good enough approximation). 

3) The effective mass density is given by Newton constant, G, times a 
quantity which may depend only on the masses of the particles involved 
(i. e. the electron mass if we study just QED effects) and the universal 
constants. 


4.4 An estimate of the dark energy density 

In the previous sections of this paper I have sketched a research program that 
in principle would allow calculating the properties of dark energy in terms of 
the fundamental quantum fields. I exclude the possibility that dark energy 
may be identified with the mean stress-energy of the vacuum fields, which 
would give an implausibly huge value. In fact I assume that a cosmological 
constant exists in Einstein equation that exactly balances that mean field 
energy, but I suppose that the cosmological constant cannot balance the 
energy of the vacuum fluctuations. Thus a calculation of the dark energy 
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properties should start from the study of the two-point correlations of the 
fluctuations of the vacuum stress-tensor. It is shown that the result would 
be finite due to metric fluctuations (see section 3). The vacuum expectation 
of the two-point correlation gives rise to a long range spacetime curvature 
as suggested in subsection 4.3, which would allow in principle getting the 
effective stress-energy tensor. It is straightforward to get a classical effective 
stress-energy tensor providing the same curvature as the two-point quantum 
correlation. Finally we should identify that effective mass density with the 
density, eq. (l2|) , of the dark energy thus getting the properties of dark energy 
from fundamental quantum fields. 

Actually a calculation along the lines of the previous subsection, involv¬ 
ing a simplified static model (i. e. a model in 3 rather than 4 dimensions), 
has been worked elsewhere [6]. The model is more realistic than the one of 
section 4.2 in the sense that it does not assume spherical symmetry. That is 
the two-point correlations of the stress-energy tensors depend on the spatial 
distance between the points, rather than on the difference between the radial 
coordinates as in the model of section 4.2. However the main conclusion is 
similar as in ours eqs. (lb111 and (1B4]) . Namely a stress-energy tensor operator 
with nil vacuum expectation gives rise to a long range spacetime curvature 
equivalent to the one produced by an effective (classical) stress-energy ten¬ 
sor.That tensor is proportional to some integral of the two-point correlation 
C{r). In particular with some plausible assumptions the said modeljB] gives 
the following effective density, p e ff, and pressure, P e ff, 



(65) 


where the numerical constant K in front of the integral is of order unity. It 
is interesting that p e ^ so obtained roughly equals the gravitational energy 
of the two-point correlations if calculated with Newtonian gravity, eq. (]5l) . If 
the vacuum fluctuations are the cause of the dark energy we should identify 
p e jf in eq. fl65|l with p DM in eq.© whence we get the following estimate for 
the QED contribution to the dark energy density 



where C (r) is identified our calculated value in eqs. (143]) and (l45|) . 

This is the QED contribution to the dark energy density, that is the contri¬ 
bution of the electromagnetic and electron-positron vacuum fields (calculated 
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to zeroth order in the perturbation). It is plausible that the contribution of 
all quantum fields would be of order eq. flSSD times a few tens. In any case the 
main obstacle for a valid estimate of the dark energy density is the difficulty 
to calculate the cut-off derived from the metric fluctuations, as discussed 
in subsection 3.2. Before solving this problem it is not worth to refine the 
calculations. 

A guess about the value of the said cut-off may be taken as follows. The 
parameter £ was introduced in section 3.1 as the inverse of the maximum 
particle’s momentum and we have argued in section 3.2 that a plausible 
estimate should be of order of some average particle mass, m, times c. Thus 
we get 

Gm 6 c 2 

Pdm = Peff ~ ^ > 

where we have written explicitly the velocity of light and Planck's constant 
for clarity. That quantity agrees with the observed value, eq.(J2]) , if we put 
m about 80 times the electron mass, which is not implausible. 


5 Excited vacuum states and dark matter 

In section 2.4 we have studied the physical vacuum state, defined as the 
ground state of the total Hamiltonian, iJ, of the interacting quantum fields, 
that is the eigenstate with the smallest eigenvalue. H should include all 
quantum fields of nature toghether with their interactions, but in this paper 
we have used the QED Hamiltonian as an example. It is a straightforward 
prediction of quantum theory that other eigenstates of H are possible in 
addition to the ground state. Amongst these states there is a particular 
class, studied in section 2.4, that I have named excited vacuum states (EVS). 
We may ask whether those states may be observed somewhere and I propose 
that they may be candidates for dark matter. 

The relevant property is that EVS consist of coupled particles possessing 
the same quantum numbers of the vacuum, that is no charges (either electric, 
baryonic, leptonic, etc.), and also no net momentum or angular momentum. 
As a consequence they would not interact via electromagnetic, weak or strong 
nuclear forces. However such states possess energy above the vacuum, see 
eq. ()38p . In section 2.5 I have assumed that the total energy of the vacuum 
state is exactely balanced by a cosmological constant, but the EVS possess 
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energy above the vacuum and therefore these states should interact gravita¬ 
tionally. 

In section 2.4 we studied EVS in a Minkowski space, that is in absence 
of graviational fields, therefore having traslational and rotational invariance. 
However it is plausible that in the presence of gravitational fields of baryonic 
matter the states are modified, suffering the atraction of large masses. A more 
detailed study would be needed in order to know the possible gravitational 
modification of the EVS. In the meantime I conjecture that the effect would 
be relevant only for extended gravitational fields like those of galaxies and 
clusters. In any case the EVS have properties not too different from those 
usually adscribed to the hypothetical constituents of dark matter, that is 
they are cold dark matter, as they possess energy but no pressure (or very 
small). A more detailed study of the possibility that dark matter consists of 
EVS will not be made in this paper, but I believe that such study is worth 
to be made. 


6 Conclusions 

The energy distribution in the quantum vacuum has been studied for the par¬ 
ticular case of quantum electrodynamics (QED). It has been shown that the 
stress-energy tensor of the combined electromagnetic and electron-positron 
fields is roughly proportional to the metric tensor, but the proportionality 
constant diverges. As a plausible solution it is proposed that there is a cosmo¬ 
logical constant in Einsein equation that exactly balances the stress-energy 
of the quantum vacuum. Arguments have been given showing that vacuum 
fluctuations give rise to a metric able to explain dark energy. Finally it has 
been suggested that dark matter may be an effect of vacuum excitations. 

7 Appendix. Two-point correlation to zeroth 
order 

The free electromagnetic field 

We shall solve the integrals in eq. (142|) , that is 
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x exp [i (k + k') • r — ek — ek'} . 

In spherical coordinates with polar angle in the direction of r we get 
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Integration of the angles <p and cp' leads to 
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where r = |r|, u = cos 9,u' = cos 9'. I 11 terms of the integrals to be defined 
below we get 
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The free electron-positron field 

We shall start from eq. (ITTj) , that is 


C D { r) = -- / d 3 p / d 3 q(E-E'f 1+ P ~ q m 

v ; 256 t r 6 ,/ P ,/ H v ; [ EE' 

x exp [i (p + q) • r —e ( p + q )}, 
that in polar coordinates becomes 
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If we approximate 
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valid for m « e -1 , we get to order 0(m 2 ) 
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The integrals that appear in the terms of order m 2 are lengthy, although 
straightforward, and they will not be reporte here. The term of zeroth order 

is 
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The integrals 

The required integrals are 


-as: 


k n dk exp (— ek ) / du cos ( kru ), 
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The results for the relevant integrals are as follows 
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